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ABSTRACT 

The distribution of global photometric, spectroscopic, structural and morphological 
iy"~) ' parameters for a well defined sample of 350 nearby galaxies has been examined. The 

, usual trends were recovered demonstrating that E/S0 galaxies are redder, more qui- 

C*") 1 escent, more centrally concentrated and possess larger Sersic indices than later type 

galaxies. 

Multivariate statistical analyses were performed to examine the distribution of all 
parameters simultaneously. The main result of these analyses was the existence of only 
two classes of galaxies, corresponding closely to early and late types. Linear discrimi- 
• nant analysis was able to reproduce the classifications of early and late types galaxies 

Oh with high success, but further refinement of galaxy types was not reproduced in the 

distribution of observed galaxy properties. A principal components analysis showed 
£^ ■ that the major variance of the parameter set corresponded to a distinction between 

^ | early and late types, highlighting the importance of the distinction. A hierarchical 

03 clustering analysis revealed only two clear natural classes within the parameter set, 

closely corresponding to early and late types. Early and late types are clearly distinct 
and the distinction is of fundamental importance. In contrast, late types from Sa to 
^ i Irr are smoothly distributed throughout the parameter space. 

' A population of galaxies classified by eye as elliptical/lenticular, and exhibiting 

concentration indices similar to early-types were found to have a significant star- 
formation activity. These galaxies are preferentially faint, suggesting they are low-mass 
systems. 
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1 INTRODUCTION 

Galaxies display a large variety of observational and phys- 
ical properties. Some properties, such as colour, emission- 
line strength and far-infrared luminosity are well correlated 
with morphological appearance su ch as bulge-to-disc ratio, 
tightness of the spiral ar ms etc. llRoberts fc Havneslll99l 
iBinnev fc M errificld 1998). These correlations are primarily 
due to the amo unt of star-forma tion, gas and dust present 
in the galaxies llKennicutdll998|) . Other properties such as 
luminosity and size, correlate well with galaxy mass, but 
have weaker correlations with morphology, since galaxies of 
different morphology have a large range in mass. 

Understanding these relationships in detail is 
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paramount to constructing a consistent picture of galaxy 
evolution. Galaxy properties are observ ed to evolve bot h as 
a function of environment and age fe.g. lBell et alll2004) . In 
order to have a clear understanding of how the observed 
evolution of each property relates to the evolution of other 
properties, it is first necessary to understand how the 
properties are related to one another. 

For example, many galaxy properties are observed to 
be dependent on the density of the galaxy's environment 
jBalogh et al.ll2004 liiogg et al]l2004ft . The morphological 
mix of galaxies changes as a function of local projected 
galaxy density, and cluster-centric radius, with the frac- 
tion of ellipticals and lenticulars increasing at higher den- 
sities, and_acOTrespondmg_decreasei fraction of spi- 
rals llDresslerl fl980l: IWhitmore et al.lll993ft . This in turn 
leads to differences in the luminosity functions of galax- 
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ies in different density regimes (e.g. ICroton et all [2005). 
The star-formation rate of galaxies declines with increas- 
ing density at densities greater than a specific threshold 
llLewis et al 1 l2002l : iBaloeh et aljEoolh . The average colour 
of galaxies becomes red der in high density environments 
l)De ProD risel^L| |2CJQ^) , again with an apparent density 
threshold (Kodama et al] l200ll) . Understanding whether 
these density relations are manifestations of the same phys- 
ical change in the galaxies, or whether several mechanisms 
are responsible, requires understanding the relations be- 
tween the morphology, colour and star- formation. 

Some galaxy properties evolve with redshift, in a man- 
ner dependent on density. At 2 « 0.5 the morphology- 
density relation shows an incre ase in the fraction of spirals 
in high density environments (Dressie r et aljll997l) . Simi- 
larly the fraction of blue galaxies in clusters increases with 
redshift iButcher fc Oemlerlll984r) . Disentangling the con- 
nected effects of galaxy age and environment is difficult, 
since transformations of galaxies due to infall into dense 
environments could produce similar observational results to 
a difference in galaxy age as a function of environment. It 
is crucial to understand the relationships between the ob- 
servational properties of galaxies at low redshift in order to 
make progress in understanding galaxy evolution. 

Many of these correl ations are alrea dy well known and 
understood. For instance iHubblel il93rJ) introduced the re- 
lation between galaxy colour and morphology. Some corre- 
lations are more recent, such as the correlation between star 
formation activity and the r\ parameter of the 2dF Galaxy 
Redsh ift Survey (2dFGRS), introduced by Madgwick et al] 
(2002), and based on a principal components analysis of 
2dFGRS spectra. 

We present here a detailed analysis of the distribution 
and correlations of the observational properties of a well- 
defined local sample of galaxies for which we have compiled 
global photometric, spectroscopic and structural parame- 
ters. We expand on previous work regarding correlations 
between galaxy properties in two major ways. Firstly, we si- 
multaneously examine correlations between a large number 
of parameters using multivariate statistical analysis tech- 
niques. Secondly, we do not limit ourselves to correlations 
with a priori assigned morphological type (discrimination) , 
though these are still investigated and form a major part of 
the work, but we also search the parameter space for natural 
clustering of galaxies (classification). 

The paper is organised as follows. After introducing the 
sample and describing the catalogue in Section [5] we de- 
scribe the characteristics of a population of blue spheroidal 
galaxies in Section [3] Following this we examine the distri- 
butions of each parameter as a function of morphological 
type in Sectional in order to gauge the overall trends dis- 
played by the galaxies. In Section|K|the parameter set is more 
thoroughly explored using a principal components analysis 
to determine combinations of parameters which account for 
a large amount of the total variance of the parameter set. 
The distribution of the first principal component separates 
early and late galaxies with good success. Section |S] further 
examines the separation of galaxies of different morphologi- 
cal type using linear combinations of parameters via a linear 
discriminant analysis. It is found that early type galaxies can 
be identified with good success, but discrimination of other 
types is less successful. In Section|7|a hierarchical clustering 



analysis is performed to look for distinct groups of galaxies 
within the parameter set, without any prior input regarding 
morphological classifications. Only two 'natural' classes of 
galaxies are found, corresponding very well with early and 
late type galaxies. The implications of the work are discussed 
in Section |H] 

A future paper will extend this work by exploring the 
results of a bulge-disc decomposition as additional input to 
the analysis. 



2 SAMPLE 

The Millennium Galaxy Catalogue 1 (MGC) is a wide, 
medium-deep, B-band imaging survey obtained with the 
Wide Field Camera on the Isaac Newton Telescope on La 
Palma (iLiske et alf2 003). It covers 37.5 square degrees with 
magnitudes in the range 13 < B < 24 mag. It is ~ 96% 
spect roscopically complete to B < 20 mag fsee lDriver et alJ 

In total MGC-BRIGHT (i.e. B < 20 mag) contains 
10095 galaxies. Although a comparatively small catalogue 
in comparison to the Two Degree Fie ld Galaxy Re dshift 
Survey (2dFGRS; ~ 250,000 galaxies. IColless et aLlliool 
and the Sloan Digita l Sky Survey (SPSS; d ata release three 
~ 370,000 galaxies. IXbazaiian etai]|2005l) . the MGC has 
several advantages that make it ideal for studying galaxy 
classification. The depth of the MGC redshift survey sur- 
passes both the 2dFGRS (6j = 19.45 mag) and the SDSS 
(r = 17.8 « B = 18.8 mag) spectroscopic limits. This in- 
crease in sensitivity allows the study of intrinsically lower 
luminosity galaxies than the larger surveys. This in turn 
allows the construction of a more representative sample of 
galaxies, and the investigation of the dependence of galaxy 
properties with absolute magnitude. 

A further important advantage is the faint surface 
brightness limit of the MGC (fiB = 26 mag arcsec -2 ). 
iDriver et al] (l2005h examine the distribution of galaxies as 
a function of magnitude and surface brightness. It is shown 
that the inclusion or omission of low surface brightness 
galaxies can have a significant effect on the resulting lu- 
minosity function. The presence of low surface brightness 
galaxies again allows construction of more representative 
catalogues. The biases resulting from a deficit of low sur- 
face brightness galaxies may be significantly different from 
biases arising from magnitude selection, since galaxies may 
be relatively luminous whilst having low surface brightness. 

Finally the quality of the images in the MGC are a 
major asset for studying galaxy classific ation. The media n 
seeing of the survey of 1.3 arcseconds (ILiske et al.l f2003). 
combined with the surface brightness limit of fiB = 26 mag 
arcsec -2 , produce very high quality images for a survey of 
this scale. The good image quality allows the structural 
properties of the galaxies to be quantified, through surface 
brightness fitting with GIM2D (Allen et al., in preparation) 
and visual inspection (Driver et al., in preparation), allowing 
this important aspect of galaxy classification to be investi- 
gated. 



1 http://www.eso.org/~jliske/mgc 
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Table 2. Number of galaxies classified into each type by SPD 
and the number of those which were agreed on by SCE. The 
percentage of agreement is also listed 





E/SO 


Sabc 


Sd/Irr 


SPD 


75 


179 


72 


SCE 


62 


140 


44 


% agreed 


83 


78 


61 



2.1 The catalogue 

We have acquired a large set of global parameters to thor- 
oughly investigate the characterstics of galaxies within the 
catalogue. The list of all parameters used in the analysis is 
given in Table Q A description, and where necessary the 
derivation, of each parameter is given below. 

2.1.1 Morphological classification 

The 350 galaxies selected for study (see Section 12. 21 have 
been visually inspected and classified into broad morpholog- 
ical classes: E/SO, Sabc and Sd/Irr. Two different techniques 
were employed by different authors to classify the galaxies, 
which makes it possible to assess the accuracy of the classifi- 
cations through comparison of the results. SPD used images 
of the galaxies printed at several surface brightness levels. 
SCE used purpose written software to display the galaxies 
alongside contours of surface brightness. 

Comparison of the classifications of galaxies by different 
authors is made in Table [5] for the 326 galaxies classified by 
SPD (with B < 19). There is » 80% agreement on E/SO 
and Sabc galaxies, but somewhat less on Sd/Irrs. The figure 
of ~ 80% may be regarded as a benchmark, against which 
to test other classification schemes. Henceforth we use the 
classification of SCE. 

2.1.2 2dFGRS ^-parameter 

The rj parameters are obtained from the final 2dF- 
GRS iColless et all l200lT) and are described in detail by 
lMadgwicl^^uTT200^) ~In brief n is derived from a prin- 
cipal components analysis of the 2dF galaxy spectra and is 
a linear combination of PCI and PC2. PCI contains infor- 
mation from emission and absorption line strength and the 
continuum in roughly equal amounts, whilst PC2 is domi- 
nated by emission and absorption line strengths. The linear 
combination is tuned to maximise the line features and is 
hence a likely indicator of the current and previous star- 
formation rate. The matching of the MGC to th e final 2dF- 
GRS catalogue is described in lCross et alJ J2004T) . No errors 
are provided for r\. 

2.1.3 Rest SDSS-DR1 Petrosian (u - r) colours 

Petrosian ugriz ma gnitudes are obtained from the SDSS- 
DR1 catalogue (see lAbazaiian et al.ll20(53) . which overlaps 
fully wit h the MGC. Details of the matching algorithm are 
given in lDriver et alJ i2005T) . as is the conversion from ob- 
served to rest-frame colour. In brief the uBgriz colours 
are used to determine the optimal spe ctral template from 
the spectral library of iPoggiantil |l998). The optimal spec- 
tra are then used to derive the appropriate K-correction 



for each galaxy through each filter. Rest-frame colours are 
derived from the rest-frame absolute magnitudes (see Sec- 
tion al,^ . The random error in the final colours are A rest- 
(u — r) ~ 0.09 mag and are dominated by the uncertainty 
in the evolutionary correction. 

2.1.4 Absolute magnitude and surface brightness 

The absolute uBgriz magnitudes and the absolute B- 
band effective surf ace brightness, are derived as stated in 
iDriver et alJ i2005f) . This includes individual K-corrections 
(as in Section 12.1.31 . a universal evolutionary correction 
(e{z) = 2.51og 10 [(l + z)-' 3 ], p = 1.5,0.75,0.75,0.5,0.25,0.13 
for uBgriz respectively), and the standard luminosity dis- 
tance corrections (assuming Qm ~ 0.3, Qa = 0.7 and Ho = 
100km s" 1 Mpc -1 ). The B-band effective surface brightness 
is derived as: \ib = m + 2.5 log 10 (27T-Re) — 101og 10 (l + z) — 
k(z) — e(z), where: m is the SExtractor measured iKronl 
(1980) magnitude, and R c is the half-light semi-major axis 
radius. The latter is derived empirically by growing the SEx- 
tractor specified ellipse until it encloses half the Kron flux 
(R° hB ) and correcting this value for the impact of the see- 
ing FWHM (T) ac cording to: R c = ^{R° bs ) 2 - 0.32r 2 (see 
IDriver et al.ll20o3 for details). The terms e(z) and k(z) are 
the evolution and K-corrections as defined above. The ran- 
dom error in the absolute B-band magnitudes is 0.06 mag 
and in the absolute effective surface brightness is 0.2 mag 
arcsec -2 . 

2.1.5 Sersic profiling 

The Sersic profiles are measured using gim2d iSimardll99ct) . 
GIm2d is a 2D image analysis code which requires the raw 
data image, a segmentation image (describing which pixels 
are associated with the object in question), and the point- 
spread function (modelled for each galaxy individually) . The 
application of GIm2d to the MGC is described in detail in 
Allen et al (in preparation). The reproducible accuracy of 
the Sersic fits have been assessed via repeat observations of 
~ 700 galaxies and the random error in log 10 n, where n is 
the Sersic index, is < 10%. The concentration and asymme- 
try pa rameters are also deriv ed via GIm2d using the defini- 
tion of lAbraham et alJ dl996>) with a = 0.1. 

2.2 Volume selected sub-sample 

For the purposes of investigating the discrimination and 
classification of galaxies into different types a sub-sample 
of the MGC was culled. This sample was selected with the 
intention of minimising bias in the properties of the galaxies 
selected whilst having a large range in magnitude, and ex- 
ploring a large parameter set of galaxy properties including 
spectral, photometric and structural information. 

There is some conflict in trying to explore a large pa- 
rameter set of galaxy properties whilst keeping selection bias 
at a minimum. Certain useful galaxy properties such as the 
2dFGRS r\ parameter (hereafter referred to simply as r\) are 
magnitude limited, since the survey from which it was de- 
rived has a limit of bj — 19.45 mag. Also derivation of pa- 
rameters such as concentration and asymmetry will likely 
be less robust for faint galaxies. Ideally a volume limited 
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Table 1. List of all parameters used in the analysis 
Parameter 
Redshift 

MGC B band absolute Kron mag., dust, K &; E-corrected (Vega mag) 

SDSS u — r absolute Petrosian colour, dust, K & E-corrected (AB mag) 

MGC B band absolute effective surface brightness, dust, K & E-corrected (Vega mag) 

2dFGR,S rj parameter 

GIM2D Sersic index 

Concentration index (inner 10%) 

Asymmetry index (background corrected) (inner 10%) 
SCE's visual classification 
SPD's visual classification 




-12 - 



0.1 0.2 0.3 0.4 



Figure 1. Absolute B band magnitude versus redshift for the 
MGC. Grey points denote galaxies for which the 2dFGRS rj pa- 
rameter is known. The sub-sample culled consists of the grey 
points within the box shown, yielding 350 galaxies. 



sub-sample would be culled containing galaxies which have 
robust measurements of all the relevant properties. However, 
due to the heterogeneous origins of some of the parameters, 
this is not possible without compromising the overall size, 
and hence the representation, of the sample. We empha- 
sise, however, that although different parameters originate 
from different surveys, any individual parameter has homo- 
geneous origins, making the distribution of each parameter 
very reliable. 

The sub-sample we selected comprised all galaxies for 
which rj has been measured within — 16 ^ Mb —22 mag 
and 0.013 ^ 2 ^ 0.05, yielding a total of 350 galaxies. The 
sample is illustrated in Figure and the catalogue is avail- 
able from the MGC web pages. Note that the sample is not 
truly volume limited since there is an apparent magnitude 
limit of bj < 19.45 mag, set by the exclusion of galaxies 
without rj, and furthermore there are surface brightness and 
size selection effects intrinsic to the MGC of /iB cff < 25.25 
mag arcsec -2 and 1.06 <half- light radius< 15 arcsecs (see 
IDriver et al"1l2005l for full details). 

We have looked for the presence of bias in the sub- 
sample, which may have been introduced by selecting only 
those galaxies for which rj is known. We have examined the 
distributions of Mb , /j,b , U — r and Sersic index for galaxies 
with and without rj within the same volume. These param- 



eters are relatively easy to interpret and should be fairly 
robust. The means and standard deviations of the distri- 
butions are given in Table [3 Also listed are the biweight 
location and scatter estimates, which may give a truer pic- 
ture of the averag e and scatter for a skewed distribution (see 
iBeers et al.lll99fj) . 

A Kolmogorov-Smirnov test was performed on each pair 
of distributions (i.e. those with and without rj) of each pa- 
rameter, to test if they are consistent with being drawn from 
the same population. The significance of each test is given in 
Table 13 where a low probability indicates that it is unlikely 
that the two sample are drawn from the same population. 

The luminosity function of galaxies in the sub-sample 
is shown in Figure [5] The solid histogram shows only those 
galaxies for which rj is known, while the dashed histogram 
shows the luminosity function for all gala xies within the se- 
lected volume. Overlaid is the best fitti ng ISch echterl (Il976l) 
function for the whole MGC JPriver et alJl2005l) . The lu- 
minosity function for the volume sele cted sample is consis- 
tent with that of IDriver et alJ <l2005l) . indicating that the 
volume is representative of the MGC as a whole. The lumi- 
nosity function for galaxies with rj has a lower normalisa- 
tion due to the fact that some galaxies have been excluded, 
but the shape is consistent, except for a deficit of galax- 
ies in the faintest bin. This indicates that the selected sub- 
sample is broadly representative of the whole MGC, except 
at the faintest magnitudes where there are too few galaxies, 
because the 2dFGRS, from which rj has been taken, only 
reaches magnitudes of 6j = 19.45 mag. This can also be seen 
in the top panel of Figure |3J which compares the distri- 
bution of galaxies for which rj is known to that of the ex- 
cluded galaxies; the excluded galaxies are generally fainter. 
A Kolmogorov-Smirnov test shows that the probability is 
extremely small that the samples with and without rj are 
equivalent. 

The middle panel of Figure |3] shows the distribution of 
rest-frame u — r colour for galaxies with and without rj. The 
distributions appear very similar. This is reflected in the 
statistics listed in Table |3] and a Kolmogorov-Smirnov test 
shows that it is reasonable to assume that the two samples 
are drawn from the same parent distribution. 

The bottom panel of Figure |H] shows the distribution 
of Sersic index for the two groups (note there are 2 galaxies 
fewer in the group without rj since the surface brightness fit- 
ting produced spurious results for these galaxies). The distri- 
butions of both groups are very similar. Again the statistics 
are given in Table [3] and a Kolmgorov-Smirnov test shows 
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Figure 2. The luminosity function of galaxies in the selected sub- 
sample. The solid histogram shows only those galaxies for which 
r\ is known, while the dashed histogram shows the luminosity 
function for all galaxies w i thin t he selected volume. Overlaid is 
the best fitting ISchechterl <1976ft function for the whole MGC 
toriver et all2005ft . 



that it is reasonable to assume that the two samples are 
drawn from the same parent distribution. 

We also checked and found no evidence for bias in the 
surface brightness of the galaxies with and without r\. This 
is at first surprising given that we are biassed against faint 
galaxies. However, whilst many of the excluded galaxies are 
faint and low surface brightness, there is a population of 
bright but low surface brightness galaxies which are included 
and hence the surface brightness distributions appear sim- 
ilar. Figure 2] shows histograms of the bivariate brightness 
distributions of the galaxies with and without r\. 

It is expected that the colours and Sersic indices of 
galaxies will approximately indicate the type of galaxy and 
underlying stellar populations. The lack of any significant 
bias in these parameters due to the exclusion of galax- 
ies without r\ suggests that the sub-sample culled provides 
a representative sample of galaxies within the magnitude 
range covered. Caution must be exercised in the interpre- 
tation of any trends in galaxy properties with magnitude 
however, since it is clear that the sample contains a deficit 
of faint galaxies. 



3 BLUE SPHEROIDS 

A histogram of the rest-frame u — r colours of the ellip- 
tical/lenticular galaxies is shown in Figure |^| The distri- 
bution is clearly bimodal, with » 60% displaying unusually 
blue colours for conventional ellipticals. This is strongly sug- 
gestive that there are two distinct populations of galaxies 
within the E/SO class. This may be due to misclassification 
of the galaxies when assigning morphologies by eye, or it 
may be revealing a fourth category of galaxies external to 
the broad classes of E/SO, Sabc and Sd/Irr. However it is 
worth noting that if due to misclassification it occurred in- 
dependently in both SCEs and SPDs samples. In order to 
ascertain the nature of the blue 'ellipticals' we have exam- 
ined the distribution of several key observational quantities 
of the galaxies. 
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Figure 3. Distributions of Mb, rest-frame u — r colour and Sersic 
index for galaxies with and without r\. The galaxies without r\ are 
preferentially faint, but have similar distributions of colour and 
Sersic index. 



The galaxies classified as elliptical/lenticular have been 
reassigned into red and blue classes based on a cut at u — r — 
2.0. Panels (a) and (b) of Figure|^]show the absolute B-band 
magnitude and surface brightness of the two populations 
versus rest-frame u — r. The blue galaxies are intrinsically 
fainter on average, with only three galaxies having Mb < 
— 18. This suggests the blue galaxies are less massive than 
the red ellipticals. 

Variations in galaxy colours can be due to variations in 
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Number Number 




Figure 4. Histograms of A/b vs. [ib for galaxies with r\ (left) and without (right). 



Table 3. Means, standard deviations, biweight average (cj,i) and scatter(st,i)i and Kolmogorov-Smirnov test significance for galaxies 
with and without r\. 

With r) (350 galaxies) Without r) (177 galaxies) K-S significance 





Mean 


(7 


Cbi 


Sbi 


Mean 


a 


Cbi 


Sbi 




M B 


-17.81 


1.10 


-17.68 


1.12 


-17.15 


1.14 


-16.67 


0.98 


2.8 x 10~ 13 


rest frame u — r 


1.74 


0.57 


1.70 


0.57 


1.67 


1.07 


1.59 


0.56 


0.24 


Sersic index 


1.70 


2.02 


1.07 


0.74 


1.58 


1.54 


1.06 


0.68 


0.88 


Ate 


22.4 


1.1 


22.4 


1.1 


22.5 


1.2 


22.6 


1.3 


0.20 



12 



10 - 



s 




1 2 3 

u-r 



Figure 5. The distribution of u — r colours for galaxies classified 
by eye as E/S0. The distribution is clearly bimodal. 

star-formation activity, metallicity, dust or age. However, it 
is unlikely that the blue population owes its colour to the 
presence of dust since this would make the galaxies redder 
than the conventional ellipticals which are known to possess 
very little dust. 

If the blue galaxies are indeed less massive than the 
red population, as indicated by their B-band luminosities, 
they will have lower metallicities, since there is a well known 
mass - metallicity relation for galaxies (see e.g. ISalzer et alJ 
I2005I for a recent analysis). Thus the mass-metallicty rela- 
tion is a plausible explanation of the variation in colour, 
since metal rich galaxies will appear redder. However, the 



bimodality cannot be due to a mass-metallicity relation for 
ellipticals /len ticulars only, since th e variation in colour is too 
large (see e.g. lStrateva et alj |2001 for the u — r colour distri- 
bution for E/S0 galaxies in the SDSS), and the luminosity 
func tion of ellipticals is wel l fitted by a Schechter function 
(e.g. lEfstathiou et al]ll98ct) and thus would not exhibit a 
bimodal colour distribution for a linear mass-metallicity re- 
lation. Thus, either the blue galaxies are not ellipticals, or 
there must either be variation in the age or star-formation 
of ellipticals as a function of mass. 

iMated lll998l) shows that all dwarf galaxies in the local 
group have an old stellar population and most also have a 
young stellar population. The relative fractions of the two 
populations and the age of the young stellar population can 
have a strong effect on the galaxy colour. Hence variation in 
star-formation activity could potentially produce a bimodal 
colour distribution as observed. 

The rj parameter of the 2d F survey is a sensitive indi- 
cator of star-formation activity ( Madg wick et alJl2002T) . The 
two classes show a strong variation in 77, panel (c) of Figure|S| 
The red galaxies occupy a narrow range of r\ with v alues typ- 
ical of massive early type galaxies (lMadgwickl2003l) . whereas 
the blue galaxies have correspondingly larger values of r\ and 
a larger scatter. Thus the red galaxies are mainly quiescent, 
but the blue galaxies display a range of star-formation ac- 
tivity. The correlation between colour and star-formation 
activity is to be expected since galaxy colour and spectra 
are intimately related, star-forming galaxies have stronger 
emission lines and display blu er colours than more quies- 
cent galaxies llKennicutdll99gl) . 

Panel (d) of Figure|S|shows the Sersic index of the galax- 
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Figure 6. The properties of galaxies classified as E/S0 as a function of colour. Blue E/S0 galaxies are preferentially faint, have higher 
r\ values indicative of ongoing star-formation and have lower Sersic indices than conventional red early-types. The concentration and 
asymmetry indices are independent of colour, suggesting that the blue ellipticals are smooth systems, in contrast to conventional late- 
types. The vertical dashed lines indicate the colour used to define the blue spheroids. The triangles mark the locations of the example 
galaxies used in Figure 171 
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ies w hen fitted with a Serisc profile only jGraham fc Driven 
120051) . The blue galaxies have values consistent with pure 
exponential systems, whereas the red galaxies have values 
in keeping with early- type galaxies. 

Although the blue galaxies have many properties in 
common with late type galaxies they have similar concentra- 
tion and asymmetry indices to normal ellipticals as shown 
in panels (e) and (f) of Figure E] suggesting that they really 
are smooth systems, and not simply misclassified late-types 
(a comparison of C with typical late types can be seen in 
panel (f) of Figure EJ. 

It seems that the blue galaxies are fundamentally differ- 
ent to both conventional ellipticals, spirals and late-types. 
Their low B-band luminosities suggests they are low mass 
systems. They have blue colours, most likely due to increased 
star-formation activity, as evidenc ed by high va lues of 77, in 
common with local dwarf galaxies ( Ma teolll998l) . We cannot 
rule out that they are misclassified late-type galaxies, e.g. 
with very low surface brightness discs, but their appearance 
and concentration suggests that they are smooth systems. 
Some examples of the appearance of these galaxies are shown 
in Figure Q along with comparison images of conventional 
red E/SOs. Hereafter these galaxies are referred to as blue 
spheroids. 



4 CORRELATIONS WITH GALAXY 
MORPHOLOGY 

The distribution of individual parameters as a function of 
morphological type has been examined. This provides an 
immediate understanding of the data set, useful for inter- 
preting the more complicated statistical analyses, and as an 
indication of which parameters are likely to form the most 
useful input into clustering, principal components and dis- 
crimination analyses. The distributions of each parameter 
are summarised below, and examples are shown in Figure [S] 
Statistics of the distributions of each parameter as a func- 
tion of galaxy type are given in Table [I] 



4.1 B band absolute magnitude 

The distribution of B band absolute magnitudes, Mb, are 
shown in panel (a) of Figure |H] The distributions for all 
classes are similar with the exception of the blue spheroids, 
which are preferentially faint, indicating that they are less 
massive systems. There is a mild trend for galaxies to be- 
come fainter going from E/SO to Sabc to Sd/Irr to blue 
spheroid, but it is less than 68% significant. 



4.2 B band surface brightness 

The distribution of B band surface brightness, fiB, is shown 
in panel (b) of Figure [S] The E/SO galaxies are generally 
bright, whilst the Sabc and Sd/Irr are ~ 1.5 mag arcsec -2 
fainter, but with a large scatter in the surface brightness of 
all types. The blue spheroid galaxies have a large spread in 
values with a peak intermediate to E/SO and Sabcd/Irrs. 



4.3 Colour 

Galaxy colours hav e been computed using u — r as 
IStrateva et alJ <l200lT) have shown that this combination can 
provide nearly optimal separation into two galaxy types for 
SDSS. Histograms of the colour distribution are shown in 
panel (c) of Figure |5] 

Recall that the distinction between E/SO and blue 
spheroid is artificial since blue spheroid galaxies are defined 
as having colours u — r ^ 2.0. There is a trend for galaxies 
to become bluer going from E/SO to Sabc to Sd/Irr to blue 
spheroid, although there is a large overlap between any two 
adjacent classes. 



4.4 2dFGRS 77 parameter 

Panel (d) of Figure |H] shows 77 as a function of type. There 
is a rather clear demarcation between E/SO galaxies and 
all the later types, and thus 7 7 shows some pro mise of a 
morphological classifier. Indeed Madgwick (2003) discusses 
in detail the ability of 77 to distinguish between early (E/SO) 
and late (spirals and irregulars) type galaxies. In that paper 
it is found that the commonly used cut of 77 = —1.4 compares 
very favourably with morphological classifications based on 
more complicated combinations of the first nine principal 
components of the spectra (77 is a linear combination of the 
first two). 

The ability of 77 to distinguish between galaxies of early 
and late morphologies is reasonably accomplished, however, 
further discrimination of different types is less robust. Panel 
(d) of Figure [5] shows a large degree of overlap between 
Sabcs, Sd/Irrs and blue spheroids. This will be shown to 
be a general feature throughout our attempts to provide a 
system of galaxy classification. It is relatively easy to sepa- 
rate E/SO galaxies from Sabcd/Irrs but further distinction 
is extremely problematic. The comparatively large spread in 
77 of the spiral gal axies is als o seen in the morphologically 
selected sample of | Madgwickl i2003T) . and in an analysis by 
iNorberg et all J2002 ) based on APM galaxies. 

iMadliwicketHd J2003T) show that 77 is much more 
tightly correlated with the rate of current star-formation 
in a galaxy (which its e lf may be correlated with morphol- 
ogy). Madgwic k et al.l i2002l) interpret 77 as a measure of 
the current star-formation rate in the galaxy, but they also 
show that 77 is correlated with morp hology for a sam ple of 
21 galaxies from the Kennicutt Atlas jKennicutJl992l) . This 
sample classifies galaxies individually into subtypes of E/SO, 
Sa, Sb and Scd. It is seen that 77 is progressively larger for 
later spiral types. Thus it could be argued that the large 
spread in our Sabc types seen in Figure |S] is a result of bin- 
ning the sub- types Sabc into a single bin. However the rather 
low 77 values of some Sd/Irrs and the overlap between the 
Sabcs and E/SOs, suggest that there is some intrinsic spread 
in the values of 77 with sub-types (this may in part be exac- 
erbated by fibre-placement issues. For example, it is possible 
that rather different 77 values would result from a fibre being 
placed on the bulge or the disc of a spiral galaxy). 



4.5 Sersic index 

The Sersic indices, n, of Sabc, Sd/Irr and blue spheroid 
galaxies fitted with Sersic profiles only (no exponential disc, 




Figure 7. Examples of galaxies classified as E/SO by eye but having much bluer colours than conventional ellipticals accompanied by 
their surface brightness profiles (top row). The bottom row shows conventional red E/SO galaxies for comparison. 
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iGraham fc Drive dl2005l) are distributed around 1, indicat- 
ing an exponential surface brightness profile. On the other 
hand E/SO galaxies have higher values of Sersic index scat- 
tered around the de Vaucouleurs value of n = 4. Panel (e) 
of Figure |H| shows the histograms. 

4.6 Concentration and Asymmetry Indices 

The concentration indices, C, are shown in panel (f) of Fig- 
urc|Hl defined as th e fraction of light w ithin 10% of the nor- 
malised radius fsee lSimard et al]l2002t) . The distribution of 
Sabcs is very similar to the distribution of Sd/Irrs showing 
an unconcentrated distribution, while the the distribution of 
blue spheroids is very similar to the distribution E/SOs and 
has higher concentration. This suggests that blue spheroids 
are not simply misclassified Sabcd/Irrs, and provides sup- 
port for the original visual classifications. 

The asymmetry indices, A, are similar for all classes of 
galaxies except Sd/Irrs which are slightly more asymmetric. 



5 PRINCIPAL COMPONENTS ANALYSIS 

The parameter set was examined in more detail by consid- 
ering groupings and correlations between several parame- 
ters simultaneously. P rincipal components analysis (PCA, 
iMurtagh fc Hecklll987l) identifies combinations of parame- 
ters which summarise the distribution of the whole set. The 
principal components are the projections onto the eigenvec- 
tors which successively contain the maximum variance, un- 
der the condition that all components are mutually orthog- 
onal. This can be thought of as finding a natural set of axes 
throughout the parameter space which minimise the spread 
of the parameters. That is, if the parameters are pictured 
as a cloud of points in space, PCA finds the axes to which 
the points are closest, where closeness is defined as the Eu- 
clidean distance between vectors. 

Thus PCA can be used to reduce the dimensionality of 
the parameter space by projecting the data onto the first few 
axes. Important parameters can be identified by examining 
the components of the eigenvectors. 

5.1 Method 

A goo d introduction to PCA is given bv IMurtagh fc Heckl 
lll987l) . The method is summarised here. 

PCA takes as its input a selection of parameters from 
the whole catalogue. This forms an N x M matrix, A, where 
M is the number of parameters and N is the number of 
galaxies. A new matrix, B, is formed by subtracting the 
mean values of each parameter from every element in the 
matrix, and normalising by the standard deviation i.e., 



Table 4. Statistics of the distributions of key parameters as a 
function of galaxy type. An idea of skewness may be estimated by 
comparing the mean and standard deviation with the biweight es- 
timators. For Gaussian distributions the values from each should 
be identical. In a heavily skewed distribution the biweight esti- 
mators are probably more useful. 
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b l3 = °'J a > = , l^i^N, l^j^M (1) 

where bij are the elements of B and Oy are the elements 
of A and 

1 N 

n j = -jyZ2 a V- ( 2 ) 

i=l 



The covariance matrix of B is then formed, 

JV 
i=l 

where Cjk are the elements of the covariance matrix C. 
Matrix C is real and symmetric, and therefore has real, pos- 
itive eigenvalues. The numerical order of these eigenvalues 
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Figure 8. Histograms of various parameters as a function of galaxy type. The median of each distribution is shown by the a dotted line. 
Panel (a) shows the distribution of Mb; (b) shows the distribution of fiB', (c) shows the distribution of u — r; (d) shows the distribution 
of r}\ (e) shows the distribution of Sersic indices; (f) shows the distribution of concentration indices. 
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yields the strength of the variance of the data in the corre- 
sponding eigenvector. 

The principal components, P, can then be found by 
projecting the eigenvectors onto the normalised, mean- 
subtracted matrix B, i.e. 

P = BC. (4) 

The order of the principal components is found from the 
numerical order of their corresponding eigenvalues, so pel 
are the elements of the A" x 1 column vector corresponding 
to the highest eigenvalue, etc. 

5.2 Results 

A PCA was carried out on various combinations of the pa- 
rameters {Mb,/xb,w — r,rj,n, C, A} in order to determine 
any particular combinations of parameters which account 
for a large amount of variance, and hence provide a useful 
description of galaxy properties. 

Many combinations of parameters were experimented 
with. In no experiment did any individual parameter dom- 
inate the principal eigenvectors: no single parameter had 
significantly more variance than any other. Considering the 
parameters as a cloud of points, they are not elongated in 
any way which is aligned to an individual parameter. Com- 
pression of the data required several parameters to be com- 
bined together, indicating that photometric (Mb, /j.b, u — r), 
spectroscopic (rj) and structural (n, C, A) parameters are 
all of similar importance in describing the distribution of 
galaxies throughout the parameter space. 

The first principal component usually projected the 
galaxies in such a way that the E/SO galaxies were sepa- 
rated from the other types of galaxies. Thus, the parameter 
space has a maximum variance in a direction that separates 
E/SO galaxies from other types. Note that PCA is not given 
any information regarding the visual morphological classifi- 
cations. The distinction between E/SO and other galaxies is 
clearly of fundamental importance in terms of their obser- 
vational properties; influencing photometric, spectroscopic 
and structural parameters. Two examples of this are shown 
in Figure fusing the subsets {Mb, fJ,B,U — r,T), n, C, A} and 
{u — r,-q,n,C}, illustrating the distinction between E/SO 
and the other galaxies. Note that the distribution of Sabcs 
is skewed with a tail overlapping with the E/SOs. Whilst 
we cannot rule out misclassification for all the Sabc galax- 
ies in the tail, some of them display very clear spiral arms. 
The other principal components showed no clear distinc- 
tion between galaxies of different types. No combination of 
principal components was found which could separate Sabc, 
Sd/Irr and blue spheroid galaxies. 

A caveat is to note that information on the errors asso- 
ciated with the parameters was not included in the analysis, 
since errors are not known for all parameters, e.g. 77. If some 
parameters have large errors associated with them they may 
skew the results. 



6 LINEAR DISCRIMINANT ANALYSIS 

Linear discriminant analysis (LDA) defines combinations of 
parameters which optimally discriminate between objects 
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Figure 9. The distribution of the first principal component, pel, 
for galaxies of different type. The top panel shows the results 
of PCA using the set {Mb, /ib> u ~ r , V, n , C, A}, and the bottom 
panel using the set {u—r, rj,n,C} .In both cases the first principal 
component clearly separates early and late type galaxies. 

of known classes. We have used LDA with various combi- 
nations of parameters to discriminate between galaxies of 
different morphology. Thus we can assess the relative merit 
of the input parameters in assigning, and hence correlation 
with, galaxy morphology. 

6.1 Method 

Consider again the N x M matrix B of equation^ contain- 
ing the normalised, mean subtracted properties of a sample 
of galaxies, where N is the number of galaxies and M is the 
number of measured properties. Assuming that the sample 
has been divided into two classes, and 1, the vectors, mo 



On the natural sub-division of galaxies 13 



and mi, describing the means of each property, and the 
covariance matrices, Vo and Vi, of each class may be cal- 
culated. 

The combined covariance matrix, W, of each class may 
then be computed, 



W = Vo+Vi. 



(5) 



Fisher's linear discriminant is defined such as to max- 
imise the distance between the samples me ans after sta ndar- 
dising by the variance within the samples llHandll98ll . page 
83). That is, we want to maximise the ratio, 



m - mi 
Wf ' 



(6) 



where f is a vector describing the set of weights which 
maximise equation |S] Thus, 



f a W (mo — mi) 



(7) 



Finally we may define a new vector containing the initial 
input data weighted so as to maximise the difference between 
class and 1, 



t = (m - m 1 ) T W _1 B. 



(8) 



Thus Fisher's LDA must first be applied to a population 
of galaxies whose types are known (in our case this is the 
whole catalogue). Thereafter an appropriate division can be 
made based on a value of t (where ti are the elements of 
vector t). The choice of t will depend on how 'pure' a sample 
is desired, setting t to give little contamination in one class, 
will give a larger contamination in the other class if the two 
classes overlap. 



6.2 Results 

The data were first standardised to zero mean and unit vari- 
ance before performing the analysis. We experimented with 
various combinations of parameters. It was relatively easy to 
separate E/SOs and the rest. Discriminating between other 
types of galaxies was always less successful. Some examples 
are given below. 

The input set {Mb, (j,b, u — r, -q, n, C, A} was used. The 
resulting discrimination is illustrated in Figure [TUl A cut at 
t(x) = 2 AO resulted in 93% of E/SOs being correctly assigned 
and 93% of others. This same cut correctly classified 99% of 
galaxies below the cut and 66% above. Thus, the separation 
of E/SO from the rest of the galaxies was relatively successful 
cf. 80% agreement when classifying by eye. 

In contrast Sabcs could only be separated with 72% 
success, Sd/Irrs with 75% success and blue spheroids with 
81% success. 

An input set consisting of purely structural parameters 
{n, C, A} was experimented with. This was less successful 
than the full parameter set, but E/SOs could still be sepa- 
rated with 83% success. 

Similarly star-forming parameters {rj, u — r} were in- 
vestigated. These two parameters alone could separate E/SO 
galaxies with 88% success, but had poorer success rates at 
separating any other types of galaxies. 




Figure 10. Example of linear discriminant analysis. The closely 
hashed histogram on the right hand side represents E/SO galaxies 
and the open hashed histogram the rest. A cut at t(x) = 2.40 
correctly classifies 92% of E/SO galaxies. 



7 CLUSTERING ANALYSIS 

So far the work presented has concentrated mainly on cor- 
relations with the morphological types assigned by eye, and 
combinations of parameters to optimise discrimination be- 
tween these a priori defined classes. It may be the case 
however, that there are natural classes of objects grouped 
within the parameter space which do not necessarily re- 
flect morphological types. A hierarchical clustering analysis 
(HCA) has been performed in order to search for such nat- 
ural groupings within the parameter set. The method and 
the results are described below. 



7.1 Method 

N galaxies, each described by M parameters, are described 
by matrix A, as in PCA. Each galaxy can be thought of as a 
point in a parameter space. The distance between each point 
and all the others is calculated yielding a set of N(N — l)/2 
distances, where the distance will be defined below. Once 
the distances are known clusters are made by agglomerating 
points within the parameter space based on their distances. 
For example, the smallest distance d^ is determined, and 
objects i and k are then agglomerated and replaced with a 
new object iUk. The set of distances are then updated using 
the new object in place of i and k. 

Several methods can be used to define distances, and 
the centres of the newly cre ated clusters. The method used 
here is that of I Ward! dl963li . which defines the distance be- 
tween two clusters as the amount by which the total distance 
from each point to the centre of its cluster would increase 
by merging the two clusters (technically, it is the sum of 
the squares of the Euclidean distances from the objects to 
the joint cluster mean minus the sum of squares from the 
objects to their individual cluster means, so the distance is 
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the amount the sum of squ ares would i ncrease if the clus- 
ters were to be merged, see lHandlll98ll) . This method was 
implemented using freely available JAVA software written by 
F. Murtagh 2 . 

7.2 Results 

Again experiments were performed with the input set 
{Mb,(J,b,u — r,ri,n,C,A}, and parameters were standard- 
ised on input. Galaxies were split into only two clearly dis- 
tinct classes. 

Classifying galaxies into two categories which gave the 
greatest separation, reveals a segregation of galaxies strongly 
dependent on u — r, r\ and the Sersic index, n. Classifica- 
tion into more than two categories did not reveal any strong 
trends with any parameter for the third class of object. Fig- 
ure shows the dependence of the two natural classes on 
u — r, r\ and n. 

The clustering of the galaxies in these natural classes 
was tested for correlations with morphology. Membership of 
group 1 or 2 was assigned as an extra input parameter into a 
LDA. No improvement (or deterioration) was made over pre- 
vious efforts to discriminate galaxies into known morpholog- 
ical types. This indicates that the classes defined from hier- 
archical clustering are no more closely correlated with mor- 
phology than the input parameters, although this could be 
due to poor classification of galaxies when assigning types by 
eye. Figure WI\ shows the distribution of each morphological 
type within the clusters. One cluster is made almost solely 
of E/SO galaxies with the majority of all other types occu- 
pying the other cluster. Although it is possible that some of 
the galaxies in class 2 of Figure [T2l were originally classified 
as Sabcd/lrr/blue spheroid as a result of misclassification, 
some of them have definite spiral morphology. 

Using the same input parameters LDA could reproduce 
the classifications of the hierarchical clustering analysis with 
excellent success of 95% (with 98% of galaxies below the 
cut, and 84% above correctly classified), showing that there 
is little overlap between the two classes in this parameter 
space. 



8 DISCUSSION 

We have examined the distributions of photometric (Mb, 
piB, u — r), spectral (rj) and structural (n, C, A) parameters 
from a well defined sample of nearby galaxies culled from 
the Millennium Galaxy Catalogue. Galaxies were classified 
into broad morphological types by visual inspection. 

8.1 Blue spheroids 

A significant fraction of galaxies classified as ellipticals were 
found to have some properties (colour, rj and n) more in 
keeping with late-type galaxies, but with concentration and 
asymmetry indices more typical of early-type galaxies. These 
galaxies are preferentially faint, suggesting that they are low 
mass galaxies. A significant fraction of these galaxies ap- 
pear to be undergoing significant star-formation. Such blue 

2 http:/ /astro. u-strasbg.fr/~fmurtagh/mda-sw/ 
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Figure 11. Galaxies classified into natural categories as a func- 
tion of rest-frame u — r, Sersic index and r]. 

spheroids have previously been reported else where, e.g. in 
studi es of faint field galaxies in HST images jDriver et alJ 
Il995l) . This population of blue spheroids may well be re- 
sponsible for the highly discrepant "early-t ype" luminosity 
function estimates fsee Ide Lapparentl 12003 for discussion). 
Early-types cuts by colour or spectral type will exclude this 
population whereas morphological cuts will not. Note also 
that the high fraction of blue spheroids in our sample is 
probably a result of the near volume-limited nature of the 
sample, the fraction would likely be much lower in a magni- 
tude limited sample. 
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Figure 12. The distribution of galaxies within each cluster per 
morphological type. 



According to hierarchical theories of structure forma- 
tion, less massive galaxies form at earlier epochs than mas- 
sive galaxies. This leads to the conclusion that the stel- 
lar populations of massive galaxies should be younger than 
those of dwarf galaxies, unless star-formation is not co- 
incident with the epoch of galaxy assembly. Observation- 
ally, however, there is a lot of evidence that less massive 
galaxies undergo star-formation at later epochs than mas- 
sive galaxies. Local dwarf galaxies are seen to be composed 
of signific ant fraction s of young stars, perhaps predomi- 
nantly so iMatecll998f) . Less massive galaxies are more likely 
to have undergone recent star-bursts than massive galax- 
ies l|Cowie et al 6; Gavazzi et al ll99dlKauffmann et alJ 
l200a) . The colour-magnitude relation of galaxies in groups 
at high redshift shows a truncation at the faint end, be- 
yond which the typical red galaxies are su pplante d by a 
blue popu lation llKod ama et all2004tlDe Lucia et al 
Similarly. iDe Propris^r^lT^OO^ shows that the increased 
fraction of blue galaxie s in clusters at high redshift (e.g. 
iButcher fe Oemlerlll984|) may be driven by low mass galax- 
ies which have recently undergone star-formation. 

The blue spheroid galaxies found in the present study 
are consistent with this picture of down-sizing in galaxy for- 
mation. The galaxies are centrally concentrated, smooth, 
and appear as normal ellipticals, suggesting they have not 
undergone any recent mergers which could trigger star- 
formation, but they possess colours and 77 values character- 
istic of star- forming galaxies. Finally it is worth noting that 
the statistical tools were unable to isolate the blue spheroid 
population from the late types. However the eyeball classifi- 
cations did. This implies that there is additional useful infor- 
mation that can be meas ured. This mi ght include smooth- 
ness (e.g., CAS ana lysis: IConselicei|2003t or low frequency 
Fourier modes (e.g.. lOdewahn^^aL 12003 for example). 



8.2 Classification and discrimination of galaxy 
types 

The over-riding feature of all the multi-variate statistical 
studies was the existence of just two clear groups of galax- 
ies (albeit with some overlap in their properties), corre- 
sponding to early types (E/SO) and the rest (Sabcd/Irr/blue 
spheroid). The physical difference being the colour, profile 
shape and the amount of star-formation activity. Further 
distinction of separate groups was always less convincing. 

The class of galaxies comprised of Sabcd/Irr/blue 
spheroids is distributed remarkably smoothly throughout 
the parameter space, and distinction between any of the 
morphological types based on their global photometric, spec- 
tral and structural properties is difficult and unconvincing. 
This is not to say there are no correlati ons of galaxy prop- 
erties with morphology. For example, iRoberts fc Havriesl 
lll994l) find correlations (notwithstanding large dispersions) 
with Hubble type and average values of total mass, mass 
fraction of neutral hydrogen, total surface density and in- 
tegrated galaxy colour, amongst others. These correlations 
are easily understood in terms of variation in star-formation 
activity as a function of morphology. However, the transi- 
tion of properties from spiral to irregular galaxies is always 
smooth, in keeping with their occupation of a single class, 
and results largely from the variation in mass along the se- 
quence. On the other hand there can be large discontinuities 
between the properties of elliptical and spiral galaxies, e.g. 
HI mass fraction. 

It is plausible then, that the astrophysical processes gov- 
erning the formation and evolution of galaxies are the same 
for all morphological types, Sa to Irr, with minor variations 
due to different masses of the galaxies, e.g. the colour differ- 
ences can be qualitatively explained by the relation between 
mass and metallicity. Ellipticals and lenticulars however, ap- 
pear to require different astrophysical processes. 

The real and fundamental separation of these two 
groups needs to be understood in the context of galaxy evo- 
lution. The existence of two distinct groups can be explained 
either in terms of evolution of galaxies or in terms of initial 
conditions. 

Perhaps the simplest explanation is that the difference 
is a result of initial conditions, with the two types of galax- 
ies having different formation mechanisms. This requires 
that ellipticals collapsed early in a dissipationless man- 
ner, and hence contain old stellar populations and a three- 
dimensional appearance. Spirals collapsed dissipatively, and 
hence formed into a disc with a reservoir of gas. The diffi- 
culty in this case lies in explaining what causes galaxies to 
collapse in different ways. Also this idea is at odds with the 
cold dark matter paradigm, which predicts that galaxies as- 
semble in a hierarch ical manner, through a series of mergers 
feaugh et aflll996j) . 

iHubbld (1936) proposed that the morphological se- 
quence of galaxies represents an evolution from the simple 
elliptical forms to the grand design of Sc galaxies. Much 
modern research proposes morphological evolution in the 
opposite dir ection, with mer ging discs forming large ellipti- 
cals (see e.g. lSchweizerl200 for a review). Reconciliation be- 
tween two clearly distinct classes of galaxies, both in terms 
of appearance and stellar population, with a transforma- 
tion between the two classes suggests that any such trans- 
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formations are rapid, or occurred at early times, in order 
to account for the relatively few galaxies with intermediate 
properties, e.g., as in Figure 

Numerical simulations show that the transformation 
from two disc galaxies to a single merged remnant can oc- 
cur in about 1.5 rotation periods of the input discs, cor- 
responding to ~400Myr for galaxies the size of the Milky 
Way dSchweizerll 20001." suggesting that morphological trans- 
formations may indeed be fast enough to preserve a bimodal 
population of galaxies. 

iBaloeh et alJ (|2004) investigate the transformation of 
galaxy colour as a function of environment, finding simi- 
larly bimodal distributions. Using lBruzual fc CharlotT ll2003l) 
synthetic stellar population models, they show that if star- 
formation is ceased abruptly then galaxy colours transform 
on the short time-scale of 750Myr. If however, star-formation 
decays with an exponential time-scale of 2Gyr, as expected 
for a process such as strangulation, they show that colours 
are transformed over a period of 4Gyr. 

Note that, for massive cluster galaxies at least, mor- 
phological transformation must take place at a slower rate 
than transformation in star-formation rate, to account for 
the existence of a tight colour-magnitude relation for lu- 
minous early-type galaxi es to high redshift iStanford et alJ 
I199SI: iHolden et aljBooi) . otherwise, at some redshift there 
would exist a population of bright blue elliptical galaxies 
which are not observed. Note also that from our visual clas- 
sifications we see very few such bright blue ellipticals, but 
we do observe apparently quiescent spirals and faint blue 
spheroids. Thus transformation of the stellar populations 
in massive cluster galaxies must be occurring faster than 
~ 400Myr if transformations are occurring at z < 1. 

Thus the bimodality of the galaxy populations could 
be due to rapid transformation of both stellar-population 
and morphology, although physically motivate d mechanisms 
such as strangulation oc cur at slower rates IlLarson et alJ 
ll980l:lBalogh et alJl200oh . 

The bimodality could be produced through slower 
transformations, if they occur mainly at high redshift. 
Ther e is indeed s ome evidence that this is the case. 
Ivan Dokkum et alJ (1999) finds a high rate of interaction be- 
tween galaxies in the z = 0.83 cluster MS1054-0 3. There is a 
close link between mergers, ultra-luminous infrared galaxies 
and quasars which suggests that the major epoch of merg- 
ing occurred c oincidentally with the peak in quasar activity, 
around z « 2 iSchweizer]l2000l) . iBundv et alJ <l2005l) present 
evidence that morphological transformations of lower mass 
galaxies may take place at later cosmic epochs, which may 
explain why the majority of blue spheroids we see are low 
mass systems. 

Note also th at the major-merger ra te, as measured in 
the MGC, is low jDe Propris et aLll2005t) . strongly suggest- 
ing that transformations as a result of mergers will not be 
very common today. Also. lLe Fevre et alJ i2000f) see a strong 
increase in the fraction of major-mergers as a function of 
redshift in the Canada-France Redshift Survey. 

These results suggest that the distinction between the 
two classes of galaxies is a result of the high redshift of the 
merging and transformation of galaxies. We note however, 
that lack of merging at low redshifts is somewhat at odds 
with predictions of CDM, in which approximately 50% of 



galaxies should hav e undergone a major merger since 2 = 1 
jMurali et afll2002l) . 

We postulate that a similar study of a high redshift 
sample of galaxies, should produce different results and a 
blurring of the two groups as the redshift of transformations 
and mergers is reached, if transformations do indeed occur 
slowly at high redshift. 

Finally we note that an an alysis of the resu lts of bulge- 
disc decomposition and CAS dConselicd l2003f) analysis of 
the galaxies are in progress and will be presented in future 
papers. 
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